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TECHNICAL PAPER 


MODELING DIGITAL CONTROL SYSTEMS WITH 
MA-PREFILTERED MEASUREMENTS 

I. INTRODUCTION 


The usual problem posed in books on digital control systems is to find a feedback controller for a 
continuous-time plant driven by a zero-order-hold with a sampled output as shown in Figure 1 (Jacquot [ 1 ] , 
p. 126). Here x(t) is an nxl state vector, n(k) is an rxl control input vector, yi(k) is an mxl output or mea- 
surement vector, F is an nxn system matrix, G is an nxr control matrix, and Q is an mxn output matrix. 
Since yi(k) = Qx(k) where k is the usual shorthand notation for time kT, yi(k) represents an instantaneous 
measure of the system at the sampling instant kT. Hence, the plant in Figure 1 will be regarded as having 
instantaneous measurements for outputs. It is well known that this system can be modeled at the sampling 
instants by the set of discrete state equations 


x(k+ 1) = A x(k) + B n(k) 

(1) 

y,(k) = C,y(k) 

(2) 


where A and B are constant matrices (Jacquot [1], p. 127). 



V|(k) 


Figure 1 . Continuous-time plant driven by a zero-order-hold with 
instantaneous measurements. 

Unfortunately, not all linear time-invariant digital control systems found in the real world have 
plants which fit the model shown in Figure 1 . For example, there exists systems in which the output, rather 
than being an instantaneous measure of the system at the sampling instants, represents an average measure 
of it over the time interval between samples. Such systems can be found in the aerospace field wherever 






startrackers and some state-of-the-art rate gyroscopes and accelerometers are used, to name a few (Polites 
[2], p. 6). In one system of this type, sets of discrete measurements are averaged to produce discrete outputs 
as shown in Figure 2. The continuous-time output z(t) is sampled every T / N seconds. Every N samples are 
averaged to generate the averaged measurement vector y A (kT), every T seconds. In Figure 2, z(t) and y A (kT) 
are both pxl vectors; C A is a pxn output matrix. Note, the plant in Figure 2 allows for the possibility of 
instantaneous measurements also. yi(kT) is an mxl instantaneous measurement vector and so Q is an mxn 
output matrix. Previously, Polites [2] developed discrete state variable formulations for the plant in Figure 
2 which have the same general form as equations (1) and (2). 



Figure 2. Continuous-time plant driven by a zero-order-hold with 
instantaneous and averaged measurements. 


Now consider the plant in Figure 3, which is a generalization of the one in Figure 2. Here, the 
continuous-time output z(t) is sampled every T / N seconds just as in Figure 2. However, now every N 
samples are multiplied by the qxp weighting matrices Hj, j = 0,...,N-1, and then summed to generate the 
qxl output vector y F (kT), every T seconds. Functionally, this is equivalent to passing the measurements 
sampled every T / N seconds through a multi-input/multi-output moving average (MA) process with coeffi- 
cient matrices H J5 j = 0,I,...,N-1 [3], The output of the MA prefilter is sampled every T seconds to 
generate v F (kT). In Figure 3, C F is a pxn output matrix and z(t) is a pxl vector. A special case exists when Hj 
= C/ N ) I, j = 0,1,. ..,N-1, where I is a pxp identity matrix. Then, the plant in Figure 3 degenerates to the 
plant in Figure 2, assuming of course that C F = C A . Another special case occurs when C F = 0. Then, the 
plant in Figure 3 degenerates to the plant in Figure 1 . 
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Figure 3. Continuous-time plant driven by a zero-order-hold with instantaneous 

and MA-prefiltered measurements. 








This brings us to a statement of the problem of this paper. Namely, what are discrete state variable 
representations for the plant with MA-prefiltered measurements shown in Figure 3? With these, the control 
system engineer can apply existing techniques to accurately design digital feedback controllers for this type 
plant. Without these, he must resort to approximate methods which could be either less accurate or else 
very time consuming if iterative. Fortunately, the latter will not be necessary, because the former is made 
possible by the formulations derived in Section III. Prior to this. Section II reviews some results to date for 
the plant with instantaneous measurements shown in Figure 1 . These will be utilized in Section III. Section 
IV presents an example which illustrates how to obtain the coefficient matrices for each formulation 
derived in Section III. Section V contains the conclusions and final comments. 


II. PRELIMINARY 


The plant in Figure 1 can be modeled at the sampling instants by the discrete state equations ( 1 ) and 
(2) where 


£ 

II 

r 

7 

3 

(3) 

A = 4*(T) , 

(4) 


and 


B = [/ 4>(\) dX] G 


(5) 


as in Reference 1 , p. 126. 4>(t) and A are nxn matrices while B is an nxr matrix. A and B can be determined 
analytically using equations (3) to (5). An alternative approach, which is also quite suitable for numerical 
computation, is as follows [4]. <}>(t) and J <|)(X) dX can be expressed in the form of matrix exponential 
series as 


<j>(t) = 



F't 1 

i! 


and 


/' <|>(X) dX = 


2 i - L - 
i = o (i+0! 


( 6 ) 


(7) 
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respectively. From equations (6) and (7), 


<f>(t) = I + F [JT <J>(\)d\] , (8) 

where I is an nxn identity matrix. Hence, / 0 T <J>(\) d\ can be determined using equation (7) with t = T and 
this result substituted into equation (8) to get 4>(T). With these results, A and B can be found using equations 
(4) and (5). 


III. STATE VARIABLE REPRESENTATIONS FOR PLANTS WITH 
MA-PREFILTERED MEASUREMENTS 


For the plant in Figure 3, 


y(t) = 4>(t - t G ) x(t 0 ) + / cj>(t — t) G u( t) dr 


(9) 


as in Jacquot [1], p. 125. Let 


v-r(k) 


.v,(k) 

.VF(k) 


(10) 


where 


.Vi(k) = Qx(k) 


( 11 ) 


and 


N-l 

y F (k) = S HjZ 
1 = 0 



( 12 ) 


Let t = kT— j(T/N) and t 0 = (k— 1)T where k = 0,1,2,... and j = 0,1,..., N— 1. Hence, equation (9) 
can be written as 
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( 13 ) 


T T kT-j(T/N) 

x(kT - j-L) = 4>CT 1)T] + J <l>(kT - j-L - T)G«(T)dT 

N N (k- 1 )T N 


Since 


z(t) = C F x(t) 


(14) 


it follows that equation (12) can be written as 


N-l T 

y F (k) = 2H j C F *(kT-j^-) 
j = 0 1N 


From equations (13) and (15), 


(15) 


N-l N-l 1 

yp(k) = [ 2 HjC F <KT-j^)]^[(k- 1)T] + 2 HjC f / 
i-0 N J -° (k- 


, kT-j(T/N) 


l)T 


^(kT-J N " T)G ^ (T)dT 


(16) 


Since 


«(t) = w[(k-l)T] , (k-l)T*£t<kT 


(17) 


equation (16) can be written as 


N-l N-l kT-j(T/N) 

y F (k) = [ 2 Hj C F 4>(T — j — )] x[(k — 1 )T] + [ 2 HjC f / <f>(kT-j-i- - T)dx]G«[(k - 1 )T] . (18) 

j = 0 N j = 0 (k - 1 )T N 

Consider the second term in the right hand side of equation (18). Using the transformation 
X = kT — j(T/N) — t in it, equation (18) becomes 


N-l N-l T-j(T/N) 

y F (k) = [ 2 Hj C F <J)(T— j— )] Jt[(k- 1)T] + [ 2 Hj C F / <t>(\) d\] G w[( k- 1)T] . (19) 

j =0 N j=0 0 
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Let 


N-l 

D + = 2 HjCp.MT-ji) (20) 

j = 0 N 


and 


N-l T-j(T/N) 

E + =[S HjC f J <|)(X)d\]G . (2D 

j = 0 0 

From equations (19) to (21), 


y F (k+ 1) = D + *(k) + E+ w(k) 


( 22 ) 


using the standard shorthand notation for kT. Letting 


-rj(k+ 1) = D+ jr(k) + E + «(k) 


(23) 


implies that 


^(k) = T|(k) 


(24) 


For the plant in Figure 3, equation (1) applies just as it does for the plant in Figure 1 . Using it and equations 
(10), (11), (23), and (24), the discrete state equations for the plant in Figure 3 can be written as 



(25) 


( 26 ) 
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A, B, D + , and E+ can be evaluated analytically using equations (3), (4), (5), (20), and (21). An alternative 
approach, which can be either analytical or numerical, is as follows. Let t = T — j(T/N), 
j = 0,1,...,N- 1 and use equation (7) to determine / T j(T/N) (j)(\)d\,j = 0,1,...,N- 1. Use these results 
in equation (8) to get <|)[T-j(T/N)], j = 0,1,...,N-1. At this point, D+ and E + can be found using 
equations (20) and (21). A and B can be evaluated using the procedure outlined in Section II. 

Notice in the formulation given by equations (25) and (26), there are (n + q) states where n is the 
number of states in x(t), and hence *(k), and q is the number of prefiltered measurements in y F (k). Hence, 
this formulation has (n + q) eigenvalues, which are the eigenvalues of A T! in equation (25). It is straightfor- 
ward to show that n of these are the eigenvalues of A and the other q are zero. 

For the special case in which Hj, j = 0,1,...,N-1 are pxp matrices such that 


Hj = 


N 


j = 0,1,...,N-1 


(27) 


where I is a pxp identity matrix, the MA-prefiltered measurements degenerate to simple averaged mea- 
surements. In this case, it follows from equations (20), (21), and (27) that 


D + — Cp A + 


where 



and 


E+ — CpB + 


where 
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These results match those previously published by Polites [2] for digital control systems with averaged 
measurements. 

For the plant in Figured, an alternative formulation to equations (25) and (26) can be derived by 
letting t = kT-j(T/N) and t G = kT in equation (9) where k = 0,1,2,... and j = 0,1,...,N-1. Hence 
equation (9) becomes 


j ( kT “j ip = 4> (-j P *0<T) + XT -0 Gu(T)dT . 

kT 1N 


(28) 


From equations (12), (14), (17), and (28), 


r N-l , . N-l kT-j(T/N) 

y F (k) = [ 2 Hj C F <}> (-jX)J x(kT) + [ 2 Hj C F f cj> (kT-jX - T ) dT]G M [(k- 

j = 0 iN j = 0 kT IN 


1)T] 


(29) 


Consider the second term in the right hand side of equation (29). Using the transformation X = kT- 
j(T/N) — t in it, equation (29) becomes 


r^ 1 / tm r N '' <.' j(T/N) l 

y F (k) = [ 2 HjC F cj)(-j jj)\ x(kT) + [ 2^HjC F / <j)(X)d\J Gw[(k- 1)T] 


i = 0 


(30) 


Let 



and 


= L 2 q HjCj 4 )(X)dxJ 


(31) 


(32) 


From equations (30) to (32), 


y F (k) = D_*(k) + E_w(k- 1) 


(33) 
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using the standard shorthand notation for kT. Letting 


T)(k+ 1) = «(k) 


( 34 ) 


implies that 


^ F (k) = D_x(k) + E_r|(k) . (35) 

From equations (1), (10), (11), (34), and (35), the alternative set of discrete state equations for the plant in 
Figure 3 is 

(36) 

(37) 


”y(k+ 1)" 


"a 0 


*(k) 


"b" 


*(k) 


= 




+ 


«( k) = [A T2 ] 


T)(k+ 1) 


0 0 _ 


_Tl(k)_ 


I 


jn(k)_ 


+ [Bj2]w(k) 



yi(k) 


"c, 0 ' 


*(k) 


.y(k) 

y T (k) = 


= 




— [Ct2J 



y F (k) 


D_ E_ 


_7)(k)_ 


T)(k)_ 


A, B, D_, and E_ can be evaluated analytically using equations (3), (4), (5), (31), and (32). An alternative 
approach, which can be either analytical or numerical, is as follows. Let t = — j(T/N), j = 0,1 ,...,N-1 
and use equation (7) to determine J^ (T/N) <}>(\)d\, j = 0,1,..., N-l. Use these results in equation (8) to get 
<M — j(T/N)] , j = 0,1,..., N-l. At this point, D_ and E_ can be found using equations (31) and (32). A and B 
can be evaluated using the procedure outlined in Section II. 

Notice in the formulation given by equations (36) and (37), there are (n + r) states where n is the 
number of states in x(t), and hence y(k), and r is the number of control inputs in u( k). Hence, this formula- 
tion has (n + r) eigenvalues, which are the eigenvalues of A T2 in equation (36). It is straightforward to show 
that n of these are the eigenvalues of A and the other r are zero. 

For the special case in which Hj, j = 0,1 ,...,N-1 , are pxp matrices as in equation (27), the MA- 
prefiltered measurements degenerate to simple averaged measurements. In this case, it follows from 
equations (27), (31), and (32) that 


D_ = C F A_ 
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where 



where 


B_ = 



-j(T/N) 

/ <t>(\)d\J G . 

0 


These results match those previously published by Polites [2] for digital control systems with averaged 
measurements. 

Now suppose the plant in Figure 3 has the prefiltered measurement vector modified as shown in 
Figure 4. From Figure 4, 


>’r(k) = y F (k)-E_«(k-l) . (38) 

Furthermore, equations (36) and (37) also apply to the plant in Figure 4. Eliminating nr|(k) and y F (k) in 
equations (36) to (38) yields the discrete state equations 


x(k+\) = 

A x(k) + B u(k) = [A T3 ] x(k) + [B T3 ] w(k) 


>’i(k) " 


’c,“ 


VT3(k) = 


= 


•j(k) = [C X3 ] x(k) . 


-' p (k)_ 


D_ 








(39) 


(40) 


Compare equations (39) and (40) for the plant in Figure 4 with equations ( 1 ) and (2) for the plant in 
Figure 1 . Notice that the plant equations (1) and (39) are the same. Only the output equations (2) and (40) 
are different. Obviously, the discrete state equations (39) and (40) for the plant in Figure 4 have n states and 
n eigenvalues, which are the eigenvalues of A. 
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Figure 4. Plant in Figure 3 with a modified prefiltered measurement vector. 


IV. AN EXAMPLE 


Consider the double integrator plant in Figure 5 where the instantaneous measurement yi(kT) 
measures x,(kT) and the MA-prefiltered measurement y F (kT) estimates x,(kT) = x 2 (kT). This example is 
rather basic in relation to the kind of problems the formulations in Section III can handle. However, it is 
chosen for pedagogical reasons to illustrate the basic procedures one uses to obtain discrete state variable 
models for plants with MA-prefiltered measurements using the formulations in Section III. 

Comparing Figures 3 and 5, 


F 


0 1 
0 0 


(41) 


G 



(42) 


C, = [1 0] 


(43) 


C F = [1 0] 


(44) 
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and 




Utilizing equations (43), (49), (50), (53), and (54), the discrete state equations for the plant in Figure 5 can 
be expressed in terms of equations (25) and (26) with 



13 


and 


B T i — 


T 2 

2 


'l T 


C T i 


I 0 0 

0 0 1 


A check will show these state equations to be observable and controllable for T > 0. 
Utilizing equations (31) and (32) with N = 2, 


D_ = H 0 C F cf>(0) + H, C F cj>( — T/2) 


(55) 


and 


r 0 -(T/2) 

E_ = L Ho C F { 4>(X.) dX + H, C F / <t>(\) d\J G 
0 0 


(56) 


Making the proper substitutions into equations (55) and (56) yields 

D_ = [0 1] (57) 

and 



( 58 ) 
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Utilizing equations (43), (49), (50), (57), and (58), the discrete state equations for the plant in Figure 5 can 
be expressed in terms of equations (36) and (37) with 


and 


1 


[A T2 ] = 


0 


0 1 0 
0 0 0 


[Bt2] — 


T 2 

2 

T 

1 


[C Tl] = 


0 


0 


0 


— T/4 


A check will show these state equations to be observable and controllable for T > 0. 

If the plant in Figure 5 is modified according to Figure 4 where E_ is given by equation (58), the 
result is shown in Figure 6. Utilizing equations (43), (49), (50), and (57), the discrete state equations for the 
plant in Figure 6 can be expressed in terms of equations (39) and (40) with 


[A t3 ] 


1 

0 


T 

1 


[B T3 ] 



T 
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and 



Again, a check will show these state equations to be observable and controllable for T > 0. One may 
recognize these to be the discrete state equations for the double integrator plant with two instantaneous 
measurements as shown in Figure 7. Hence in Figure 6, x 2 (kT) is being observed by sampling Xi(t) every 
T/2 seconds, prefiltering the samples with a backward difference equation, and modifying the result as 
shown in Figure 6. 






V. CONCLUSIONS 


Two different discrete state variable formulations were derived for the continuous-time plant driven 
by a zero order hold with a combination of instantaneous and MA-prefiltered measurements. This plant is 
shown in Figure 3 and the discrete state variable formulations are presented in equations (25), (26), (36), 
and (37). The first formulation has (n + q) states where n is the number of states in y(t), and hence ,y(k), and 
q is the number of prefiltered measurements in v F (k). The second formulation has (n + r) states where r is 
the number of control inputs in w( k). If the prefiltered measurement vector is modified as shown in Figure 4, 
a third formulation can be derived as shown in equations (39) and (40). This one has n states and the same 
plant equation as equation ( 1 ) for the plant in Figure 1 . 

With regard to which formulation is the best to use for modeling a given plant with MA-prefiltered 
measurements, the following advice is offered. Choose the one which yields the least number of states in 
the resulting discrete state equations. In general, this is the formulation which has n states, equations (39) 
and (40). However, this formulation does require modifying the prefiltered measurement vector as shown 
in Figure 4. If this is either impossible or else undesirable, for whatever reason, then choose the formula- 
tion, from the remaining two, which yields the fewest number of states. Hence, if the plant in question has 
fewer prefiltered measurements than control inputs (i.e., q < r), choose the formulation which has (n + q) 
states, equations (25) and (26). If the plant has fewer control inputs than pre filtered measurements (i.e. , r < 
q), choose the formulation which has (n + r) states, equations (36) and (37). If the number of prefiltered 
measurements equals the number of control inputs (i.e., q = r), the choice is arbitrary. By following these 
guidelines, the resulting set of discrete state equations is more likely to be observable and controllable, plus 
there will be fewer states and, hence, fewer eigenvalues to contend with in designing the feedback 
controller. 

Having the formulations derived in this paper, the control system engineer can apply standard avail- 
able design techniques to accurately design a digital feedback controller for any plant which fits the model 
shown in Figure 3. He does not have to resort to approximate methods which may be either less accurate or 
else time consuming, if iterative procedures are required. 
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